Optical Bloch oscillations in periodic structures with metamaterials 
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We predict that optical Bloch oscillations can be observed in layered structures with left-handed 
metamaterials and zero average refractive index where the layer thickness varies linearly across the 
structure. We demonstrate a new type of the Bloch oscillations associated with coupled surface waves 
excited at the interfaces between the layers with left-handed material and conventional dielectric. 



Electron oscillations in the presence of a constant elec- 
tric field were predicted by Bloch in 1928 [lj. Such Bloch 
oscillations become possible due to beating of the local- 
ized eigenmodes of the structure corresponding to the 
equidistant eigenstates of the spectrum known as the 
Wannier-Stark ladder [2J. Experimental verification of 
the theory was impossible at that time, since dephasing 
time of electrons in crystals is shorter than the period of 
the electron Bloch oscillation. Later, electron Bloch oscil- 
lations were observed in semiconductor superlattices [5] 
for which the period was reduced due to a small mini- 
band width in the artificial structure. 

Dephasing processes for electromagnetic waves are neg- 
ligible making the observation of the optical Bloch oscil- 
lations in photonic systems much easier. The first ex- 
perimental observation of optical Bloch oscillations was 
reported in Ref. [I] for linearly chirped Bragg gratings. 
Later, several studies reported the observation of optical 
Bloch oscillations in various structures [5j [fH [3 19] . 

Recent experimental realization of left-handed ma- 
terials |10j has opened up many unique opportunities 
to explore novel effects in the structures with nega- 
tive refractive index. In this Letter we study, for the 
first time to our knowledge, optical Bloch oscillations 
in one-dimensional layered structures containing alter- 
nating layers of left-handed and conventional dielectric 
slabs. We choose the material parameters in such a way 
that the average refractive index n across pair of the 
neighboring layers vanishes, thus fulfilling the condition 
for the existence of a novel type of the specific zero-fi 
bandgap [TTJ [12]. We change the layer thickness lin- 
early in the structure and observe an optical analogue 
of the Wannier-Stark ladder in the eigenmode spectrum, 
and the corresponding Bloch oscillations in the resonant 
transmission bands. We reveal that in such structures 
the Bloch oscillations can be observed in three different 
regimes. Compared to the photonic Bloch oscillations 
in conventional dielectric structures, the metamaterial 
structures can support a novel type of the Bloch oscil- 
lations associated with coupling of surface waves at the 
interfaces between left-handed and dielectric layers. 

We study a one-dimensional layered structure shown 
schematically in Fig. 1 , where the slabs with the width bi 
are made of metamaterial being separated by a dielectric 
slab with the width a*. Variation of the refractive index 
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FIG. 1: Schematic of linearly chirped one-dimensional pho- 
tonic crystal with alternating layers of left-handed metama- 
terial and dielectric. 



in the i-th pair of layers can be described as follows: 
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where ni and n r are the refractive indices of metamaterial 
and dielectric, respectively. 

We consider TE-polarized waves with the electric field 
having one component E = [E X ,Q, 0), and waves propa- 
gating in the plane (y, z). In this case, the field distribu- 
tion can be described by the Hclmholtz equation: 



A 2 E x (y,z)+n 2 (z)E x (y,z)- 



1 dfj,(z) dE x (y,z) 
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dz 



dz 



= 0, (2) 



where A2 is the two-dimensional Laplacian, and the co- 
ordinates are normalized to c/lu. Firstly we consider 
periodic structure. Electric field in an infinite one- 
dimensional periodic structure can be represented as a 
superposition of Bloch eigenmodes [13], with the electric 
field envelopes U(z + A) = U(z), where A is the struc- 
ture period. The dispersion relation for the Bloch waves 
is found by the transfer matrix method |13j . 
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FIG. 2: (Color online) Bandgap diagram for the TE-polarized 
waves. Black and white areas correspond to gaps and bands, 
respectively. Two spectra of the excited Bloch oscillations are 
shown on the left. The inset shows a magnified part of the 
spectrum. 



where Kb is the Bloch wavenumber, k z 



~Fy n fr ~ an d k y is the normalized propagation con- 
stant along the y axis. According to this relation an 
infinite stack containing metamaterials exhibits a non- 
resonant gap for n = A -1 j n{z)dz — 0, where n is an 
average refractive index of the structure, i.e. n r a = \nib\. 
This condition is easy to fulfil for negative refractive in- 
dex materials. 

The bandgap diagram of the layered structure is shown 
in Fig. 2 for the parameter plane (A, k y ). Here we assume 
that dielectric is vacuum, e r = /i r = 1, and that it is two 
times thicker than the second layer, a/b = 2. We choose 
the parameters of the left-handed media as follows: e r — 
—5 and u r = —0.8. This set of parameters allows surface 
waves to exist at the interfaces between metamaterial 
and vacuum [IS]. As follows from Fig. 2, for the zero 
n structure the bandgap spectrum differs substantially 
from the case of conventional periodic structures made 
of conventional dielectrics [13] • Stack with the average 
zero refractive index possesses a complete gap with the 
transmission resonances [Til H2] when the optical path of 
the wave in either layer of the period coincides with a half 
of the wavelength in the corresponding medium. Thus for 
the normal incidence (k y = 0) transmission is observed 
only when n r a = nib = Ttm, where m is integer. For the 
slabs of equal thickness the regions of the transmission 
resonances in (A, k y ) plane degenerate into infinitely thin 
lines. 

We study the propagation of electromagnetic waves in 
such a layered structure with zero average refractive in- 
dex in each pair of layers, when the thickness of layers 
is chirped linearly, i.e. A q = Aq + qSA, where integer q 
numbers the layers. We are looking for localized solu- 
tions in the structure, and for numerical simulations we 
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FIG. 3: (Color online) Field distribution in the case of surface- 
wave-assisted Bloch oscillations. The Wannier-Stark ladder 
appears for the propagation constants centered around k y o = 



2.47, normalized period is L y 



820. 



consider a finite stack of layers with perfect metal bound- 
ary conditions, E(z — 0) = E(z = L) = 0, where L is the 
total length of the structure. We assume the Gaussian 
field distribution in the plane y — across the layers, and 
in order to find the electromagnetic field distribution in 
the whole stack we look for its eigenmodes by solving the 
Helmholtz equation Then we decompose the initial 
field distribution using the basis of eigenmodes and find 
the solution in the whole structure. 

To find the eigenmodes of Eq. ^ we employ the fol- 
lowing discretization scheme |14j : 
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where z rn = mh are the mesh points with the discretiza- 
tion step h and E x (y, z) — U (z)e~J kyV . Such a discretiza- 
tion scheme provides an algorithm convergence |14j , and 
it avoids excitation of spurious modes in the structure. 

In metamaterials the energy flow /[Ex H]dz can be 
negative, i.e. the energy can propagate in the opposite 
direction to the propagation constant k y [TS]. Conse- 
quently, we determine the direction of the energy flow 
of each eigenmode and choose the sign of the propaga- 
tion constant such that the energy flows in the positive 
y-direction. Decomposition of the initial condition in the 
plane y = in the eigenmode basis is made using the 
least squares method. 

To find propagation constants (values of k y ) which lead 
to the Bloch oscillations, we analyze the spectrum of 
eigenvalues of this layered structure. The Bloch oscil- 
lations are expected to appear where the spectrum of 
eigenmodes is equidistant. Practically for all gradients 
of a linear ramp we observe several sets of equidistant 
states. The equidistant eigenvalues of k y correspond to 
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FIG. 4: (Color online) Field distribution for the case of guided 
waves. The Wannier-Stark ladder appears for the propagation 
constants centered around k y o = 1.34, period is L y = 100. 

a spatial optical equivalent of the Wannier-Stark ladder 
which is associated with the Bloch oscillations. 

Spectrum of k y can be divided into three different re- 
gions. First, when k y < n r < electromagnetic waves 
propagate in both left- and right-handed materials. In 
the second region, n r < k y < \m\, waves propagate 
in metamaterial only being evanescent in the vacuum 
layers. In this regime, our structure can be considered 
as an array of coupled left-handed waveguides. When 
k y > \ni > \n r , only surface waves may propagate along 
the interfaces separating different materials. 

We find that the Bloch oscillations can be observed 
in all three regimes of the wave propagation when the 
corresponding set of equidistant propagation constants is 
excited. We consider a stack containing 36 pairs of meta- 
material and dielectric slabs and the normalized period A 
varying from 3.7 to 6. First, we excite the eigenstates cor- 
responding to the regime of surface waves with the cen- 
ter of the spectrum at k y0 — 2.47. Figure 3 presents the 
intensity distribution for the electric field which shows 
clearly spatially periodic oscillations of the beam posi- 
tion in the structure. The corresponding spectrum of 
eigenstates is shown on the left side of Fig. 2. We note 
that the beam reconstructs its shape after each period of 
oscillations. The field is highly confined to the interfaces 
between metamaterial and vacuum, demonstrating that 
such Bloch oscillations exist due to interaction of surface 
waves in the structure. The distance between Wannier- 
Stark eigenstates Ak y defines the period of oscillations, 
L y = 2ir/Ak y . For this case, we find L y = 820, and this 
agrees well with Fig. 3. 

Bloch oscillations of the beam with the spectrum cor- 



responding to the coupled waveguide regime, n r < k y < 
\ni\, are shown in Fig. 4. The equidistant spectrum 
of eigenstates corresponding to the Wannier-Stark lad- 
der is also shown in Fig. 2 (top, left). We notice that 
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FIG. 5: (Color online) Field distribution for the case when 
waves propagate in both type of materials. The Wannier- 
Stark ladder appears with average propagation constant 
k y o — 0.8, period of oscillations is L y — 210. 



oscillations are strongly anharmonic, but they are still 
periodic with the period defined well by the relation 
L y = 2ir/Ak y , which is less than the period of Bloch 
oscillations associated with surface waves. 

The regime of Bloch oscillations corresponding to the 
waves propagating in both media can be found in a 
different structure with wider transmission resonance. 
We analyse that the structure consisting of 36 periods 
and where the normalized period varies linearly from 
2.5 to 7.5 (corresponding to the period change gradi- 
ent 5 A = 0.14). Ratio of the layer thicknesses in each 
period is the same as in the previous calculations. We 
choose e = —3.6 and fi = —1.11, preserving the zero av- 
erage refractive index of the structure. The calculated 
field distribution in this case is shown in Fig. 5, and the 
center of the equidistant spectrum appears at fc y o ~ 0.8. 

In conclusion, we have studied the propagation of 
electromagnetic waves in layered structures with left- 
handed metamaterials, and have demonstrated that lin- 
early chirped structures with zero average refractive 
index can support novel types of Bloch oscillations. 
We have demonstrated that the excitation spectra are 
equidistant, manifesting a similarity with the optical 
Wannier-Stark ladder. Using numerical simulations, we 
have demonstrated three different types of the Bloch os- 
cillations, and we have revealed the existence of unusual 
oscillations associated with coupled surface waves. 
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